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Abstract 

A procedure for tailoring a blended wing-body 
configuration to reduce its computed wave drag is 
described. The method utilizes an iterative algo- 
rithm within the framework of first-order linear the- 
ory. Four computed examples are included. In each 
case, the zero- lift wave drag was reduced without an 
increase in drag due to lift. 

Introduction 

In the initial rough-cut state of supersonic con- 
figuration design, the requirement for low wave drag 
is a prime consideration. Often, at this stage in the 
design process, linear codes are used to compute the 
lift and the drag due to lift (ref. 1) as well as the 
zero-lift drag (ref. 2). 

Reference 2 includes a design option for tailoring 
the fuselage of a configuration to obtain low zero- 
lift wave drag. However, this procedure has several 
weaknesses. It is applicable only if the configura- 
tion has a fuselage that is defined by circular cross 
sections. The procedure is based on the zero-order- 
accuracy Eminton-Lord theory (ref. 3), by which the 
computed drag is independent of Mach number. This 
independence is inconsistent with the first-order the- 
ory (ref. 4) that is used in the analysis phase of the 
program to compute drag from the equivalent area 
distributions obtained from Mach plane slices. Fi- 
nally, the procedure does not account for the volume 
that is required to fill the gap that occurs between the 
fuselage and the wing when these two components are 
input separately in the wave-drag geometry format. 

This paper describes a method for tailoring a 
blended wing-body configuration for low wave drag. 
The method utilizes Mach plane slice area distribu- 
tions and so is consistent with the wave-drag analysis. 
Furthermore, since for a blended wing-body the fuse- 
lage geometry is not input as a separate component 
to the wave-drag program, there is no requirement 
to account for the fuselage-wing gap. The method 
designs for low drag at zero lift in accordance with 
the first-order analysis procedures of reference 2, and 
then the methods of reference 1 are used to ensure 
that this tailoring does not increase the drag at the 
design lift coefficient. 


Ci lift coefficient 

Cd,wO zero-lift drag coefficient 

D w zero-lift wave drag 

E(x ) error distribution 

L length of equivalent body 

M Mach number 

MSAD Mach slice area distribution 

q dynamic pressure 

S synthesized distribution obtained by 

averaging individual equivalent area 
distributions 

S r cross-sectional area of body of 

revolution 

S@ equivalent area — area of Mach 

plane slice at polar angle 6 through 
configuration, projected onto 
Y Z-plane 

T thickness along centerline section 

V volume; volume parameter 

X = x/L 

x,y,z Cartesian coordinates 

A D W { wave-drag component corresponding 

to 0i set of Mach plane slices 

77, £ dummy integration variables 

0 angle between F-axis and projec- 

tion of Mach plane slice normal vec- 
tor onto Y Z - plane 

Analysis 

Linear Wave-Drag Equations 

Reference 4 shows that the linear theory ap- 
proximation to the wave drag of a general, slender, 
nonlifting configuration is 


Symbols 


B 

base area of body of revolution 

B b 

equivalent base area corresponding 
to displacement effect of jet exhaust 

Cd,vjI 

drag coefficient due to lift at design 
lift coefficient 



where is the wave drag associated with the 

0 set of Mach plane slices. These are the Mach 
planes whose surface normal vectors project onto the 
FZ-plane at angle 0 to the F-axis. (See fig. 1.) 



Equation (1) is approximated by the sum of a finite 
number of Mach slice drag components: 

1 n 

D u , (2) 

i= 1 

The increments A D w j are defined by the formula 

AD w,i = -^ J o J o [ S 8i (0 S 0 i (v) Ink - v\ dt, dr] 

(3) 

where Sgj is the equivalent area associated with the 
Mach plane slices at angle The primes indicate 
differentiation with respect to x. The equivalent area 
is the area of the projection of the Mach plane slice 
through the body onto the TZ-plane (fig. 1). 

The Design Problem 


sometimes a net increase. This method succeeded in 
reducing the drag only for poor initial designs (those 
having very high drag); consequently, it was rejected. 

Detailed Procedure 

The method that was adopted utilizes a more 
global approach. The procedure is described here, 
and a flow diagram is given in figure 2. 

First, the individual equivalent area distributions 
are computed and averaged. The result is a single 
synthesized distribution having a value of zero at the 
nose (x — 0) and terminating in a value B which is 
the equivalent area associated with the displacement 
effect associated with the jet plume and the wake. 

Next, this synthesized distribution is compared 
with the minimum drag area distribution for a body 
of revolution having a positive base area B = 

(ref. 7): 


The design problem for low wave drag is to find 
a feasible configuration geometry for which the sec- 
ond derivative distributions in equation (3) yield a 
low value of D w in equation (2). It may be instruc- 
tive to compare this problem with the design prob- 
lem for low sonic boom (ref. 5), which specifies a 
target equivalent area distribution for only one set 
of Mach slices (that at 0 = —90°). With the plan- 
form and the camber distribution considered to be 
fixed, a straightforward iterative procedure can be 
used to tailor the thickness distribution to obtain the 
required equivalent area distribution (ref. 6). In con- 
trast, the wave-drag design problem involves multiple 
equivalent area distributions. Tailoring the thickness 
distributions for any one of these could have a neg- 
ative impact on the other distributions. This effect 
was demonstrated in an initial abortive attempt at 
synthesizing a design procedure. For this attempt, 
16 Mach plane sliced area distributions were used. 
Because of symmetry only nine distributions needed 
to be calculated. Each distribution was written out, 
together with the wave drag associated with that 
distribution. The distribution associated with the 
largest value of wave drag was then subjected to a 
smoothing procedure to generate a highly smoothed 
version of the original distribution. The smoothed 
distribution was then used as the target distribution 
for tailoring the thicknesses of the wing-fuselage sec- 
tions to reduce the wave drag for this particular Mach 
slice distribution. This procedure always reduced the 
wave drag for that distribution and, to a lesser ex- 
tent, it usually reduced the wave drag associated with 
the adjacent distributions. For the remaining distri- 
butions the drag increments varied between positive 
and negative so that the total drag increment was 


Sr (x) 


SV - B 


(l - X 2 ) m + -Xy/l - X 2 + - cos " 1 X 

7 7 r 7 r 

( 4 ) 


where V is the volume of the body of revolution; but 
w r hen equation (4) is used as a target distribution for 
an aircraft geometry, V is treated as a design param- 
eter that controls the configuration volume, which 
is generally slightly less than V. An option in the 
computer code described in reference 2 permits one 
to compute the actual overall volume of a configu- 
ration. The value of V to be used in equation (4) 
to obtain the required configuration volume is ulti- 
mately determined by iteration. 


The third step in the procedure is to alter the 
thickness distribution of the wing-body along Mach 
cone slices. Thus, for each value of x for which the 
synthesized distribution is computed, a double Mach 
cone is constructed with its axis aligned with the 
flight direction axis and its vertex at the given value 
of x. For this j:, the synthesized distribution is com- 
pared with the target distribution. If the synthesized 
distribution is larger, the thicknesses are reduced 
proportionately where the Mach cone intersects the 
wing- body. Similarly they are increased where the 
synthesized distribution underestimates the target 
distribution. These thickness variations can be ap- 
plied along the forward part of the Mach cone, the 
rearward part, or both. 


This process is iterated until an approximation 
to the target distribution is obtained. The target 
distribution will not be obtained precisely because 
the thickness changes are applied vertically (in the 
z-direction) and not in the Mach slice direction. 
Consequently, attempts to refine a design beyond 
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a reasonable point tend to cause “wiggles” in the 
distribution area. These wiggles result in large values 
in the second derivatives in equation (4) and thereby 
increase the drag. Only three or four iterations 
are required to obtain a smooth approximation to 
the target distribution (fig. 3). One determines the 
actual number of iterations required by examining 
the drag computed by equation (2) in each analysis 
iteration, stopping the iteration when it begins to 
increase, and backing up one iteration. 

The final step in the procedure is to satisfy the 
volume constraint. The volume of the redesigned 
configuration is compared with the required volume, 
and the value of V in equation (4) is altered propor- 
tionately. Then the entire procedure is repeated to 
obtain the required volume within a specified error 
bound. 

We evaluate the design by comparing the zero-lift 
wave drag of the configuration with that of the orig- 
inal configuration using 64 Mach slice distributions 
in the analysis. Also, the induced drag due to lift at 
the required lift coefficient is compared with that of 
the original configuration. 

There is little mathematical basis for the above 
procedure. If the drag increment computed by equa- 
tion (3) were a linear functional of the correspond- 
ing area distributions, then the synthesized average 
area distribution could be used to compute an aver- 
age drag. As it is, this average area distribution can- 
not be used to compute any meaningful drag value 
except for a body of revolution. 

Intuitively, however, the procedure does provide a 
means of obtaining a kind of global improvement in 
the volume distribution so that the sum of the indi- 
vidual drag increments is decreased. The procedure 
has succeeded in reducing the computed wave drag 
for all cases attempted thus far, even when consider- 
able effort has been expended to generate a low-drag 
initial configuration. Several of these sample cases 
are described in the following section. 

Sample Cases 

The following examples were generated strictly 
for the purpose of illustrating the effectiveness of the 
procedure for reducing the wave drag. No attempt 
was made to satisfy the many constraints associated 
with factors such as engine sizing or distribution of 
weight and volume that would be required for a prac- 
tical configuration designed for a specific mission. On 
the other hand, some effort was made to generate 
initial configurations with relatively low wave drag, 
so that a further reduction of the drag would be a 
significant accomplishment. Although the thickness 


was tailored to reduce the drag, the planform was 
not changed, and the overall volume, the lift, and the 
base area associated with the jet exhaust were held 
essentially constant. All the configurations included 
a vertical fin and four generic engine nacelles but 
no horizontal tail or canard (fig. 4(a)). For the sec- 
ond example, the fuselage cross-section shapes were 
nearly circular, but for the other cases, the fuselage 
cross sections were more nearly elliptical (fig. 4(b)). 

Figure 5 shows the planform for the first example: 
a configuration designed for flight at M = 1.6. 
The overall planform area is 11 686 ft 2 , the fuselage 
length is 295 ft, and the wing span is 155 ft. The 
wing leading edge is continuously curved, with a tip 
sweep angle of approximately 60°. The thickness 
distribution along the centerline is shown in figure 6, 
along with the revised thickness distribution. The 
results of the thickness tailoring are as follows: 


Parameter 

Initial 

Revised 

C'D.wi) 

0.00286 

0.00132 

Cd.wI 

.00583 

.00580 

c L 

.10033 

.10001 

v; ft 3 

40 467 

40 395 


The zero-lift wave drag is considerably reduced 
without significant variation in the other parameters. 

Figure 7 shows the planform for the second ex- 
ample: a Mach 2.0 configuration. The planform area 
is 12160 ft 2 , the fuselage length is 295 ft, and the 
wing span is 156 ft. The leading edge has a straight 
segment having a sweep angle of 66°. The centerline 
thickness distributions are shown in figure 8, and the 
design results are as follows: 


Parameter 

Initial 

Revised 

C'D.wG 

0.00269 

0.00245 

O D,wl 

.00733 

.00712 

c L 

.09831 

.09816 

v, ft 3 

41826 

41814 


In this case, the reduction in Cjy >w 0 IS on ^y a bout 

10 percent. This is augmented by a slight reduction 

in Cd,wI’ 

Figure 9 shows the planform for example 3, an- 
other Mach 2.0 configuration. The planform area is 

11 975 ft 2 , the fuselage length is 295 ft, and the wing 
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span is 155 ft. The primary leading-edge sweep an- 
gle is 67°. The centerline thickness distributions are 
shown in figure 10, and the design results are as fol- 
lows: 


Parameter 

Initial 

Revised 

Cd,w 0 

0.00236 

0.00175 

C D.wl 

.00667 

.00644 

C L 

.09945 

.10005 

V, ft 3 

40127 

40176 


In this case a significant decrease in Cq w q is 
obtained, together with a slight decrease in Cq w i . 

The fourth example is a Mach 2.4 configuration 
whose planform is shown in figure 11. The planform 
area is 12 278 ft 2 , the fuselage length is 295 ft, and 
the span is 151 ft. The primary leading-edge sweep 
angle is 69.1°. The centerline thickness distributions 
are shown in figure 12, and the design results are as 
follows: 


Parameter 

Initial 

Revised 

Cd,w 0 

0.00182 

0.00154 

Cd.wI 

.00610 

.00603 

Cl 

.08620 

.08577 

V, ft :} 

40 212 

40235 


Even though the Cq w q value for the original con- 
figuration is already quite low, it is further reduced 
in the design process. 

For the sonic-boom-type designs (first, third, and 
fourth examples), a thinning occurs toward the aft 
portion of the section, but a compensating thickening 
occurs forward of this region, so that a nearly con- 
stant volume is maintained. This type of thickness 
redistribution occurs because of the far-aft maximum 
span position of the wing. This revised thickness dis- 
tribution may pose additional structural and space- 
allocation problems. In the second example, how- 
ever, which has the wing in a more forward location, 


the revision results in a more equitable thickness 
revision. 

Concluding Remarks 

A procedure for tailoring a blended wing-body 
configuration to reduce its wave drag has been de- 
scribed. The procedure utilizes an iterative loop 
within an analysis code that computes the zero-lift 
wave drag from multiple Mach plane slice area dis- 
tributions. In the design process the geometry is con- 
strained by holding essentially constant the planform 
shape, the overall volume, the lift, and the base area 
associated with the jet exhaust. Four sample cases 
were described. In each case the zero-lift wave drag 
was reduced without an increase in drag due to lift. 

NASA Langley Research Center 
Hampton, VA 23681-0001 
August 3, 1992 
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Figure 2. Flowchart of design procedure. 

















(a) Initial. 



(b) Revised. 



(c) Comparison of centerline thickness distributions. 

Figure 6. Thickness distributions along centerline for Mach 1.6 configuration. 











Figure 9. Planforin of Mach 2.0 sonic-boom configuration. 
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(a) Initial. 



(b) Revised. 



(c) Comparison of centerline thickness distributions. 

Figure 10. Thickness distributions along centerline for Mach 2.0 sonic-boom configuration 










(b) Revised. 



(c) Comparison of centerline thickness distributions. 

Figure 12. Thickness distributions along centerline for Mach 2.4 configuration. 
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